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The recent theoretical discovery of families of travelling; wave solutions in pipe flow 
l|Faisst fc Eckharddl2003tlWedin fc Kerswellll2004HHof et alEoO^ at Reynolds numbers 
lower than the transitional range naturally raises the question of their relevance to the 
turbulent transition process. Here a series of numerical experiments are conducted in 
which we look for the spatial signature of these travelling waves in transitionary flows. 
Working within a periodic pipe of 5 D (diameters) length, we find that travelling waves 
with low wall shear stresses (lower branch solutions) are on a surface which separates 
initial conditions which uneventfully relaminarise and those which lead to a turbulent 
evolution. Evidence for recurrent travelling wave visits is found in both 5 D and 10 D 
long periodic pipes but only for those travelling waves with low-to-intermediate wall 
shear stress and for less than about 10% of the time in turbulent flow. Given this, 
it seems unlikely that the mean turbulent properties such as wall shear stress can be 
predicted as an expansion over the travelling waves in which their individual properties 
are appropriately weighted. Rather, further dynamical structures such as periodic orbits 
need to be isolated and included in any such expansion. 



1. Introduction 

Wall-bounded shear flows are of tremendous practical importance yet their transition 
to turbulence is still poorly understood. Typically, the laminar flow solution is linearly 
stable (e.g. plane Couette flow, pipe flow) or if linearly unstable only well beyond the 
regime where transition occurs (e.g. channel flow). As a result, transition is an abrupt 
process triggered by a perturbation of sufficient amplitude. Gcnerically, this could be ex- 
pected to lead to an intermediate state of reduced symmetry but in fact the flow always 
immediately becomes temporally and spatially complicated. A recent new direction in 
rationalising this phenomenon revolves around identifying alternative solutions (beyond 
the laminar state) to the governing Navier-Stokes equations. In the past few years, such 
solut ions in the form of steady states or travelling waves have been found in plane Couette 
flow llNagat^lQQOHClever fc Busselll992l Il997t IWaleffelll998h. channel flow llltano fc Tohl 

mil 



l200lUWaleffeil200ll l2003l). an autonomous wall flow Jimenez fc Simensll200l|) and most 



recently pipe flow ijFaisst fc EckhardtlEool IWedin fc KersweUlbOO^T lnvariablv. these 
solutions are saddle points in phase space. The idea is that through their stable and 
unstable manifolds and interconnections between them together with the appearance of 
periodic orbits, phase space becomes sufficiently complicated to support 'turbulent' tra- 
jectories. 
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Gathering supporting evidence for this scenario is in its infancy especially for spatially- 
extended systems but progress is being made. Eckhardt and coworkers have been central 

in app l ying such dynam i cal sy s tems ideas to transit i on to t urbulence fe.g.lSchmiegel fc Eckhardtl 
4l997|) : lEckhardt et all (|2002j) : lFaisst fc Eckhardtl l|2004fk lEckhardt et all (|2007|) ). Their~ 
basic approach has been to focus on the statistics of many transition events rather than 
specific examples and so far have been reasonably successful in rationalising the proper- 
ties of reduced numerical m odels with known dynamical systems structures in phase space 
(see lEckhardt et all (|2007|) for a nice review). Complementary work has concentrated on 
establishing connections between specific flow behaviour a nd underlying nonli near solu- 
tions present. In channel flow, for example, Itano and Toh l|ltano fc Tohll200l|) interpret 
wall turbulent 'bursting' events with flow along the unstable manifold of a travelling wave 
solution. They also managed to isolate a periodic-looking solution on the basin boundary 
of the turbulence by c ontinually adjustin g a tr ajectory such that i t neither relaminarised 
or became turbulent l|Toh fc Itanoll2003|) (se elSkufca et all l|2006|) for an equivalent cal- 
culation in a model system). Jimenez et al. l|^U^enez^^zll2005|i studied both channel 
flow and plane Couctte flow in an effort to relate near-wall turbulent events to the large 
number of known nonlinear solutions. They concluded that the turbulence stayed close 
to the upper branch travelling waves as far as comparing simple statistics of the flow 
field such as maximum (over space) wall-normal and streamwise components were con- 
cerned. Other work has focussed upon identifying isolated periodic solutions directly 
from numerically-integrated turbulent t rajectories using a Newt on-Raphson technique. 
In both the case of pl ane Couette flow llKawahara fc Kidall200ll) and highly-symmetric 
forced box turbulence Ijvan Veen et aZJl200(Tl) , the authors claim to find one periodic orbit 
which seems to share the same mean properties as the turbulent attractor. 



In pipe flow, the only work so far aiming to establish th e physical relevance o f the 
recently-discovered travelling waves has been experimental l|Hof et allfoOQd. 12005^ . By 
analysing the flow structure acros s turbulent pipe flow (both of 'puff' and 'slug' type - see 
IWvgnanski fc Champagne] 1 197.^ . good correspondence was found, at least occasionally, 
with the outer symmetrically-arranged ring of fast 'streaks' (streamwise velocity anoma- 
lies) which is one of the dominant features of the travelling waves. The match is less clear, 
however, with regards to the complementary slow streaks centred arou nd the pipe axis as 
well as with the smaller cross-stream velocities (e.g. figs 2E and 2F of lHof et all 1 20041) 1. 
The purpose of this paper is to build on this work by carrying out a detailed quantitative 
study which can explore how closely the travelling waves are reproduced or 'visited' in 
phase space and the frequency of such visits using direct numerical simulations. If it 
emerges that turbulent pipe flow can be understood as an effectively random switching 
between the neighbourhoods of these travelling waves then an appropriately weighted 
expansion across the 'active' travelling waves visited may provide a useful predictor of 
the turbulent flow properties. This presumes that some version of periodic orbit theory 
develope d in low dimensional dynamical systems (e.g. ijCvitanovid Il988t lArtuso et all 
may carry over to this very high (formally infinite) dimensional setting. 

The structure of the paper is as follows. §2 begins by briefly describing the numerical 
method used to solve the Navier-Stokes equations before discussing reasonable measures 
chosen to quantify if, and how well, the flow approaches a travelling wave (TW) solu- 
tion. There is a certain amount of arbitrariness in this choice because the TWs are fully 
nonlinear solutions not obviously orthogonal under any inner product. Hence some exper- 
imentation has been necessary before a final choice on the exact 'correlation' functions to 
evaluate made. Given also the fact that the TWs are parameterised continuously by their 
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axial wavelength (albeit over a finite range), it has been convenient to impose a strict 
periodicity in the pipe to discrctisc the TWs which can exist in the system. As a result, 
a periodic pipe of length 5D has been used for the majority of the results. Even then, 37 
TWs of two-, three- and four-fold rotational symmetry about the axis can be found at a 
Reynolds number of 2400. These are briefly described in §3 together with their stability. 
In §4, we show numerical evidence that some of the TWs are visited but not all and for 
only part of the time. In §5, the statistical frequency of these visits is quantified by exam- 
ining the correlation data from across a number of runs. Finally, a discussion follows in §6. 



2. Formulation 

2.1. Numerics 

The Navier-Stokes equation and solenoidal condition for the flow of an incompressible, 
Newtonian fluid along a circular, straight pipe under the action of an imposed pressure 
gradient are 

d t u + u.Vu+ -Vp = vV 2 u, V.u = 0, (2.1) 
P 

where v is the kinematic viscosity, p the pressure and p the constant density. Non- 
dimensionalising the system using U the mean axial speed and the pipe radius D /2 where 
D is the diameter gives rise to the Reynolds number Re := UD/v. A constant mass flow 
flow rate - or equivalently Re - is maintained along the pipe at all times. A numerical 
solution for the primitive variables (velocity and pressure) was developed in cylindrical 
coordinates (s, 9, z) using finite differences in the radial direction (s) and Fourier modes 
fo r z and 9. The time stepping was performed using the third-order Runge-Kutta scheme 
of lNikitinl <|200fTl) . 

Hereafter a quoted numerical resolution of (N, M, K) corresponds to N + 1 equally- 
spaced radial points (i.e. a grid step of 1/N where < s < 1) and Fourier expansions in 
9 and z of wavenumbers —M/2, ...,M/2 and (-K/2, K/2)ir / L respectively where L 
diameters is the nominal length over which periodicity is imposed. As discussed below, 
the main choice of pipe length was L = 5: for this geometry a coarse grid was (25,32,30), 
an intermediate grid (50,48,40) and a fine grid (50,60,60): equivalent grids in longer pipes 
were used (e.g. a fine 10D-grid was (50,60,120)). 

In addition to checks of specific components of the code using analytic test solutions, 
a series of calculations was performed using as initial conditions a TW solution plus a 
perturbation in the form of the leading (unstable) eigenfunction. Good agreement was 
obtained for the growth of the disturbance and that predicted from the eigenvalue. Im- 
portantly, the code was also cross- validate d with another time-stepping code based on 
velocity potentials ijWillis fc Kerswellll2006l) . 

Commonly in studies of this type the grid is stretched in the radial direction as the 
highest resolution is required near the wall. However, here there was little difference be- 
tween the results obtained with uniform and non-uniform grids. This is not surprising 
given the nature of the flow at the Reynolds number considered in this study, as can be 
seen in Figure below. 
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2.2. Travelling Waves 

The TWs so far identified l|Faisst fc Eckhardtl2003tlWedin fe KerswellfcOO^ are arranged 
into symmetry classes of m-fold rotational symmetry about the axis and then continu- 
ously parametrised by their axial wav enumber across a finite range. TWs of 1-fold through 
to 6-fold symmetry have been found l|Kerswel]ll2005() but only 2-, 3- and 4- fold TWs are 
currently known to exist below Re = 2485. Within each symmetry class, the TWs ap- 
pear through saddle node bifurcations so that close to the saddle node point there is a 
well-defined upper and lower branch solution for a given wavenumber. For higher Re, 
the solution surfaces typically kink and fold back on themsel ves so that multiple pairs of 
branches can coexist at the same wavenumber (e.g. figure 10 o flWedin fe Kerswei]|l)2004|) '). 

Imposing a pipe periodicity immediately reduces the continuum of TWs present down 
to a discrete number which can fit into the pipe. This is a crucial simplification which 
means the matching procedure adopted below can monitor all the TWs available to the 
flow. The main choice of a 5D long pipe was a compromise between the need to keep the 
number of TWs to a manageable size (helped by a shorter pipe) and the need to have a 
dynamical system which could support turbulent behaviour at a value of Re where the 
TWs ar e fully resolvable (helped by a longer pipe). T his pipe length has also been studied 
before feggels et q^Jll994t iFaisst fc Eckhardtll2004|) and sustained turbulence predicted 
for Re > 2250 (IFaisst fc Eckhardtll2004h . This short le ngth, of course, precludes captur- 
ing turbulent spatiotemporal features such as 'puffs' ijWvgnanski fc Champagnelll973(l 
which typically extend over 20-D but does allow an examination of 'temporal' turbulence 
which, when triggered, fills the whole pipe. 

Figures Q]and[5] show the result of tracing out all the 2-, 3- and 4-fold solution branches 
in the friction fa ctor-axial wavelength plane at Re — 2400 where the friction factor 
JSchlichtinglll968l) is defined as 

1 dp I U 2 , . 

A := f- / (2.2) 

pdz/ 2D y ' 

with dp 1 1 dz bei ng the mean pr essure gradient. The curves are similar but less contorted 
at Re = 2000 l)KerswelJl2005|) and Re = 2200 (not shown). The vertical dotted lines 
drawn at axial wavenumbers a = 0.625n where n = 1,2, ...5 indicate the TWs which fit 
in the pipe which is actually taken to be 7r/0.625 = 5.0265 diameters long. The m-fold 
symmetry class, the letter label and wavenumber are used to identify the TWs in what 
follows. For example, the TW with wavenumber 1.25 and lowest friction factor in figure^ 
is the 3&_1.25 TW. For this pipe geometry and Re, there are 37 TWs (6 2-fold, 22 3-fold 
and 9 4-fold rotationally symmetric TWs) which can be numerically resolved and used 
to match against the flow. 



2.3. Matching 

As fully nonlinear solutions, the TWs do not possess any simple orthogonality with 
respect to an inner product. Therefore establishing when a directly numerically simulated 
flow, udns, approaches a chosen TW velocity field, utw, is not a straightforward case 
of projection. Given this, a number of ad-hoc 'correlation' functions were developed and 
tested to measure how close the flow comes momentarily to a TW. The basic approach was 
to construct an integral based on the velocity fields u^ns an d W-tw over one wavelength 
27r/a of the TW near the middle of the pipe. The correlation functions found most useful 
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Figure 1. Solution branch for 3- fold rotationally-symmetric travelling waves plotted on a fric- 
tion factor versus axial wavenumber plot at Re = 2400. The blue dotted line represents the 
lower bound given by the Hagen-Poiseuille solution (A = 64/ Re) and the upper blue dashed 
line corresponds to the Re = 2400 value of the log-law parametrisation of experimental data 
= 2.0 log(Re m y/A) — 0.8 (see Schlichting 1968 equation (20.30)). The solution branch is only 
shown as far as it is assured resolved (hence the loose ends: the main mapping resolution was 
(8, 30, 6) in the truncation nomenclature of Wcdin & Kerswell 2004). The dotted vertical lines 
indicate the wavenumbers (a = 0.625n in units of 2/D, n = 1,2,3,4,5) which fit into a pipe 
of length 7r/0.625D long. The letters are used to label each allowable TW together with the 
wavenumber. 



and thereby adopted were a normalised inner product based on the total velocity fields 



Itot(t) = max 
once the mean profile of the TW 



(2.3) 



Wtw{s) ■= -j-2 / / w TW (s,0,z) dOdz (2.4) 



4tt 2 Jo 

had been subtracted from both velocity fields 



v_djvs := u DNS (s, 6*o + 0, z Q + z) - w T w{s)z 
vtw ■= utba(s, 6, z) - wtw(s)z 
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Figure 2. The equivalent of figure but for the 2-fold (left, blue) and 4-fold (right, black) ro- 
tationally-symmetric travelling waves. Typical truncations used to resolve this solutions were 
(9, 25, 7) for the 2-fold case and (6, 40, 5) for the 4-fold case. 



and an inner product using the cross-stream velocity components 



Iuv(t) 



max 

#0,^0 



DNS' U TW 



DNS' U DNS 



TW U TW 



Here u 1 - = (u,v,0) is the cross-stream velocity part of u, and 



,ui,u 2 ) := — j 



2-rv/a r 2ir 



Ui • U2 sds d6 dz 



(2.5) 



(2.6) 



is the (usual) inner product. 

The second measure I uv is necessary as the TW streamwise (streak) velocities are 
typically an order of magnitude larger than the cross-stream velocities and hence their 
matching contribution tends to dominate Itot, to the extent that there is little difference 
between I to t and I w (the equivalent of I uv but using the streamwise velocity component 
only), even though the mean flow has been subtracted. A correlation using the full ve- 
locity fields vdns and vtw is not a sensible measure as it is dominated by the mean 
streamwise flow to the extent that all correlations involving the axial velocity would be 
large. 
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By design, Itot and I uv can only take values in the interval [—1, 1] with a value 1 
indicating a perfect match. The phase optimisation over 9o and zq (carried out by sys- 
tematically evaluating all the options over the [0, 2n) x [— 7r/a, ir/a) grid) in practice 
ensured that the correlations were never very negative typically lying in the interval 
[—0.2,0.2]. Experience indicated that there is evidence for a TW visit if I to t and I uv 
reach values of 0.5 and above (although more on this below). 

A number of other measures were tried. An inner product using the streamwise vortic- 
ity, which provides a single measure of the cross stream flow, was found to shadow I uv , 
although at a lower level. The symmetry assumed for the TW's (see 12.71 below 1 ) ensures 
that the cross stream velocity is zero along the axis of the pipe. This would not be ex- 
pected to occur over an extended length of the pipe in fully turbulent flow. Hence an inner 
product over a subset of the domain excluding the central portion could be suitable. This 
was investigated, and it was found that restricting the domain to i < s < 1 produced 
somewhat higher higher correlations, particularly for I uv , but, again, a similar pattern of 
behaviour. Hence, these other measures could be used to produce essentially the same re- 
sults by adjusting the level of correlation that would be regarded as giving a good match. 

The matching was performed by maximising the value of over all possible values 
of 9q and Zq. Itot was then calculated for the same orientation. I to t is not suitable as 
the primary measure as it is so heavily dominated by the streamwise component that 
the cross-stream structure of the flow would, in effect, be discounted when choosing the 
"best" match. 

Some of the TW's are highly correlated (values greater then 0.95 have been observed), 
so that, at a specific time, using these measures, there can be a good match of the flow 
to more than one TW. In cases such as this, the value of the mean wall shear stress and 
the perturbation kinetic energy can be used to select the most appropriate match. 

2.4. Travelling Wave Stability 

An obvious way to start the DNS runs is to use a TW together with some small pertur- 
bation as an initial condition. If this perturbation is unstructured, for example, by relying 
on numerical discretization errors, the flow takes a long time to exit the neighbourhood 
of the TW and wastes cpu time. A better strategy is to find the unstable eigendirections 
of the TW and to use the most unstable cigenfunction with a small amplitude as the 
perturbation. Four 'lower' branch TWs - 2b_1.25, 3a_2.5,3h_2.5 and 4b_3.125 - and four 
'upper' branch TWs - 2a_1.25, 3b_3.125, 3j_2.5 and 4c_3.125 - were selected as starting 
TWs. The distinction between upper or lower branch solutions can be ambiguous when 
the solution surface is as convoluted as in figure ^ Here we consider that TWs with high 
friction factors are upper branch solutions and those with low friction factor are lower 
branch solutions. The 8 choices made represent extreme and therefore unambiguous ex- 
amples under this categorisation. 

The stability properties of the 8 chosen TWs are listed in Table 1. The travelling waves 
all possess the shift-&-reflect symmetry S 

S : (s,9,z) -> (s,-0,z + ir/a), S : (u,v,w,p) -> (u, -v,w,p). (2.7) 

so permitted linear disturbances can either be partitioned into those symmetric or an- 
tisymmetric with respect to S. When checked the TWs were invariably stable to an- 
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Branch TW Number of Unstable Largest Growth Rate Resolution 

Eigenvalues (in units of U/D) (M, Af, C) 



Lower 2b_1.25 lr 1.1 x KT 1 (9,25,7) 

3a_2.5 2r 2.1 x KT 1 (8,30,6) 

3h_2.5 3r+2c 1.7 x 1CT 1 (8,30,6) 

4b_3.125 lr+2c 2.6 x lO" 1 (6,40,5) 

Upper 2a_1.25 4c 6.6 x 10~ 2 (9,25,8) 

3b_3.125 2c 3 x 10" 3 (8,30,5) 

3j.2.5 6c 1.7 x 10" 1 (8,30,6) 

4c_3.125 6c 3.7 x 10" 1 (6,40,5) 



Table 1. The stability properties of typical upper and lower travelling waves at Re = 2400 : r 
and c indicate the number of real and complex eigenvalues respectively. The resolution (M,N , C) 
is the same for the travelling waves and the stability calculation and indicates the azimuthal, 
radial and axial resolution respectively: see Wedin & Kerswell (2004) for details. The unstable 
eigenvalues all correspond to disturbances possessing the same shift-&-reflect symmetry as the 
travelling wave. 



tisymmetric disturbances so Table 1 concentrates exclusively on the situation in the 
iS-symmctric subspacc. The number of unstable directions is strikingly small given the 
large degrees of freedom involved (e.g. 0(15,000)) and the size of the growth rates - 
O(0.1 U/D) indicate incrtial instabilities. Figure |31 shows the structure of the most un- 
stable eigenfunctions for the lower branch TWs 26_1.25 and 3a_2.5. The 3a_2.5 TW gives 
a particularly clear example of how an unstable eigenfunction is concentrated in the re- 
gions of maximum shear in the TW streak velocity. 



3. Results 

A series of runs were performed by taking as initial conditions each of the selected 8 
upper and lower TWs perturbed by a small amount of their most unstable eigenfunc- 
tion. Since this perturbation can be added or subtracted, 16 runs were in fact done. 
This protocol highlighted a fundamental difference between the upper and lower branch 
TWs. For all the lower branch TWs tested, starting the run in one sense along the TWs 
most unstable manifold invariably led to an uneventful gradual rclaminarisation, whereas 
starting in the other sense always produced a turbulent evolution. Both signs of pertur- 
bation, in contrast, produced a turbulent trajectory for the upper branch TWs. This 
implies that the 4 lower branch TWs (and by implication other lower branch TWs) and 
their stable manifolds are part of a boundary dividing regions of phase space which lead 
to the two different types of behaviour. Their unstable manifolds are then either directed 
towards the laminar or turbulent states normal to this bounding surface. The terminol- 
ogy 'basin of attraction' and 'separatrix' has purposely been avoided here as the second 
initial observation when doing these runs at Re = 2400 with a 5D pipe is that the turbu- 
lence is only transitory. In other words, the laminar state is still the global attractor at 
Re = 2400 in a 5D periodic pipe although the flow can experience a long but ultimately 
finite turbulent episode. Not surprisingly, the length of these turbulent transients seems 
to depend sensitively on the exact numerical resolution used. This is particularly true 
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Figure 3. The travelling waves 2b_1.25 (top left) and 3a_2.5 (bottom left) and their most 
unstable eigenfunctions (top right for 2b_1.25 and bottom right for 3a_2.5). The arrows indicate 
the cross-stream velocities (larger arrows corresponding to larger speeds) for all 4 plots. The 
shading represents the axial velocity for the eigenfunctions whereas for the travelling waves, the 
axial velocity differential away from the laminar flow corresponding to the same mass flux is 
plotted (dark most negative - slow streaks - and light most positive - fast streaks). The same 
contour levels are used throughout with the eigenfunctions renormalised so that their largest 
axial velocity is set to the largest absolute contour level (the shading outside the pipe indicates 
0: contours levels range from —0.417(7 to 0.266(7 in 8 steps). 



with regards to the azimuthal and axial resolution - reducing this too far can produce 
what looks to be sustained turbulence (e.g. turbulence remains transient when reducing 
the resolution from the working resolution (50, 60, 60) to (25, 30, 30) or (50, 24, 24) but 
lo oks to be sustained at (50 , 16, 16) over a time greater than 3,000D/J7: the resolution 
of lFaisst fe Eckhardtl l)2004[) which predicts sustained turbulence at Re = 2250 is some- 
where in between these last two choices). The issue of exactly when (or indeed if!) pipe 
flow turbulence becomes sustained is an area of much current in terest at the moment 
l|Peixinho fc Mullinll2006t iHof et q/.lE006l IWillis fc KersweUll2006[) . 



The two runs started around the lower branch TW 46_3.125 illustrate the general 
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behaviour well. Figure 01 plots the kinetic energy (per unit mass) 



KE:=— / / i(u-u iam ) 2 sdsdOdz (3.1) 
J Jo J a - 

of the flow beyond the laminar state versus the mean wall shear stress 
i r 2L r 27T n w 

T:= 2^LDj I PU ~d7 _ 1 Mdz l=- 1 s PU 2 A = -±ReA(2 P U 2 /Re)} (3.2) 

for the flow evolution together with all 37 of the TWs present. Strictly when comparing 
with values for a TW, the kinetic energy and mean wall shear should be calculated for 
the section of the pipe containing the best match. However, for a pipe of this length, 
flow structures tend to persist over the full length of the pipe, particularly in the near 
wall region containing the streaks. The wall shear stress depends only on the strcamwise 
velocity, which also provides the dominant component of the kinetic energy perturbation. 
As a result, there is usually little difference in these quantities for the full pipe and a 
section corresponding to one of the travelling waves. Larger differences were observed in 
longer pipes. 



For one sign of the eigenfunction perturbation, the flow tamely relaminarises while for 
the other it executes a long turbulent transient. During this latter evolution, there is 
evidence of close visits to at least 5 TWs. The first (labelled '1' in figure^} is to 4/_3.125 
and occurs during the early stages as the flow trajectory moves away from 46J3.125: see 
figure El This quality of the match is extremely high suggesting that there may be a het- 
eroclinic connection between the two TWs. The fact that even at t = there is already 
a considerable correlation with 4/_3.125 indicates that 46_3.125 is structurally similar to 
4/_3.125, an observation which is also true for some other groupings of TWs within the 
same rotational symmetry class. Figure also shows the correlation signal for 4c_3.125 
which has a very close visit after 120 D/U (labelled '4'). Taken together, the correlation 
functions for the other 4-fold TWs indicate that the flow retains its 4-fold symmetry 
until about 130 D/U (this was verified by examining the transient solution), whereupon 
it switches to a predominantly 3-fold symmetry. Figure shows this switch-over well via 
the correlation function I tot for 3a_3.125. This TW is the best candidate for a close visit 
(labelled '5' at t ~ 250 D/U) over all the 3-fold symmetric TWs. Figure shows how 
the contribution to the correlation functions vary over the matching wavelength at point 
'5'. The total correlation I to t is uniformly high but is low indicating that there is 
probably a good match in the streak structure but not in the cross-stream velocities. 
This is confirmed in figure where there is clear evidence of 3 equally-spaced fast streaks 
around the outside of the pipe like 3a_3.125 which are fairly streamwise independent. 
Surprisingly, there is little similarity between the inner slow streak structures. The axial 
vorticity is a convenient but particularly discriminating way of probing the cross-stream 
velocity fields as it involves taking derivatives. Given the low value of I uv , the poor 
comparison is expected however there is evidence that the DNS flow and the TW have 
roughly the same wavelength of variation. 

During the dominantly 4-fold symmetric phase, there is also evidence for visits to 
26_0.625 and 2a_1.25 (labelled '2' and '3') shown in figure |H1 A further way to charac- 
terise the extent of all these visits is to compare the kinetic energy and wall shear stress 
of the flow across the length of the pipe used for matching with the values associated 
with the TW. This data is collected together in figure 01 The relative closeness of these 
points to the corresponding values for the TW (e.g. '5' and 3a_3.125) is further support- 
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Figure 4. The surplus kinetic energy per unit mass, in units of U 2 versus wall shear stress r 
in units of 2pU 2 /Re for udns starting at TW 46.3.125 (marked as a o with ax). The solid line 
indicates the turbulent evolution for one sign of the eigenvalue perturbation and the thick dashed 
line traces out the uneventful relaminarisation for the other (respectively 46_3.125(+/— )). The 
laminar state is represented by the point (—4,0). All the TWs present are also plotted: A for 
2-fold TWs, + for 3-fold TWs and o for 4-fold TWs. Filled in symbols indicate TWs which 
appear to be visited by udns and the numbered dots indicate the temporal points of closest 
approach. In chronological order: 1 - 4/_3.125 (lower black filled circle); 2 - 26.0.625 (blue filled 
triangle furthest from 3); 3 - 2a A. 25 (closer blue filled triangle); 4 - 4c_3.125 (upper black filled 
circle) and 5 - 3a_3.125 (circled red +). 

ing evidence of a visit. 

Runs started at upper branch TWs with either sign of unstable eigenfunction lead to 
turbulent-looking trajectories. Figure ITU1 shows one of these runs starting with 36.3.125. 
Just as in the case of 46.3.125, the initial rotational symmetry of the flow lingers for a 
substantial time. During this phase, there is clear evidence of close visits to 36.2.5, 3c_2.5 
and 3j_2.5 (see figure ITT|) . The first visit is particularly significant as the structural over- 
lap between the initial TW (36.3.125) and the visited TW (36.2.5) is low: by point '1', 
both I tot and I uv for 36.2.5 have risen from w 0.4 to near 0.8 in just over 50D/U. 
Figure H~51 which compares the initial condition, the DNS flow at point '1' and the TW 
36.2.5, confirms that the flow makes a significant adjustment to match with the new TW. 

After a time of about 125 D/U - see Figure [T3*l- the flow switches to being more 2-fold 
symmetric. The flow then successively visits 2a_1.25 and 26.1.875 before relaminarising 
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Figure 5. hot and I uv (inset) as a function of time for TWs 4/_3.125 (red solid line) and 
4e_3.125 (blue dashed line) starting at TW 46_3.125. Dots label the times of likely closest visits 
to each TW (times '1' and '4' coincide with those in figure QJ. 

at a time 350 D/U. Point 4 gives a very good comparison between the the DNS flow 
and the TW 2a_1.25 because the correlation functions I to t and I uv arc simultaneously 
large within the comparison wavelength: see Figure ITU at z m 1.5. The velocity plots - 
see Figure 1151 - are very similar even up to the same vortex in the top left quadrant of 
the cross-sectional snapshots. 

In Figure ITUI point 5 is flagged as a match to TW 26.1.875 as this TW has the largest 
correlation to the flow at this point, with I to t = 0.63 and I uv = 0.61. However, this is 
not the only TW that satisfies the criterion for a good match. Specifically, TW 2a_1.25 
has Itot — 0.58 and I uv = 0.49. Also, 2a_1.25 is closer to point 5 than 26_1.875 in Figure 
1101 That both of these TW's show a high correlation with the flow is not surprising as 
the TW's also show a high correlation (matching 26_1.875 to 2a_1.25 gives I to t = 0.53 
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Figure 6. hot and I uv (inset) as a function of time for TW 3a_3.125 starting at TW 4b_3.125. 
The dot labels the time of closest visit (corresponding to '5' in figure 0. 

and I uv = 0.62). Figure IT7)1 shows a comparison between the flow at point 5 and the two 
travelling waves. In general, Figure El suggests that there is better agreement between 
TW 2a_1.25 and the flow than 26_1.875, despite the fact that it has a lower correlation 
as measured by I to t and I uv . That this can occur is not unexpected since hot and /„„ 
are measures of shape and not amplitude. Figure ^] and the position of the TW's on 
1101 suggests there is a better agreement in amplitude between 2a_1.25 and the flow than 
with 26.1.875. 

No significant correlations were found with any 4-fold TWs as the correlation func- 
tions I tot < 0.3 and I uv < 0.45 throughout the evolution. This was a persistent finding 
in all the other runs where a 2-fold or 3-fold symmetric TW was used to initiate the 
flow. Even when the flow was started by a 4-fold symmetric TW, correlations with other 
4-fold symmetric TWs would only be significant in the initial phase of the flow evolution 
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Figure 7. The left plot shows the correlations I to t (upper line) and I uv (lower line) over one 
wavelength of TW 3a_3.125 at point 5 in figure|S] The right plot shows the azimuthal distribution 
of the wall shear stress in units of 2pU /Re at the axial position of maximum I to t + Iuv near 
z = 1.5. The upper lines give the azimuthal stress and the lower lines the axial stress (minus 
the laminar value of —4). The regular lines with 3- fold symmetry are for the TW and the more 
irregular ones from the DNS values. 

where the seeding symmetry is still present. In contrast, certain 2-fold and 3-fold TWs 
were consistently visited regardless of how the flow was initiated. The fact that the 4- 
fold symmetric TWs have a higher wall-shear-stress-to-kinctic-encrgy ratio than the flow 
seems to adopt is indicated by figures 0] and EH and equivalent plots for other starting 
TWs (not shown) . The flow never seems to visit the part of phase space where the 4- fold 
TWs are unless specifically inserted whereas some of the 2-fold and 3-fold symmetric 
TWs are in the active part of phase space populated by turbulent trajectories. 

Information about which TWs are visited in each of the 12 turbulent runs is sum- 
marised in Tables El and 01 The criterion used to indicate a visit is Itot > A and 
Itot + Iuv > 2A with A = 0.5 and A = 0.6 for a closer visit. There is a strong corre- 
lation between the rotational symmetry of the starting TW and the waves subsequently 
visited even when the initial transient - defined as the period when I to t for the initial 
TW decreases (typically of O(30 D/U) in duration) - is not considered. For instance, runs 
1,2,5 and 10 - see TablcEl- show no evidence of visits to TWs with different rotational 
symmetry. However, across this suite of runs, all the TWs are visited at least once except 
for the 3 TWs 3<x_1.25, 36_1.25 and 3c_1.25. Tabic G3 shows the more discerning results 
of only considering times when the flow leaves the rotational symmetry class of the ini- 
tial TW. Examples of when this occurs have already been discussed in run 4 (at about 
130 D/U: see FigureEJl and in run 7 (at about 125 D/U; See Figure 113) . This makes it 
particularly clear that the 4-fold rotationally symmetric TWs are never visited except 
when the run is specifically started near one of these TWs (runs 4, 11 and 12) whereas 
2-fold and 3-fold symmetric TWs are visited regardless of the starting symmetry (e.g. 
runs 3, 6, 7, 9). 



3.1. Frequency of Visits 

Given the evidence above that the TWs are visited, the next question is how frequently. 
To answer this, 'visit' statistics were compiled across a number of different runs which 
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Figure 8. Comparison plots of the DNS flow (left column) and the TW 3a_3.125 (right column) 
at point 5 in figure The top row shows the velocity fields at the streamwise position of 
maximum I to t + /«» shown in figure [7] (CHECK) The shading represents the axial velocity 
perturbation from laminar flow with contours from -0.55 (dark) to 0.5 (light) for the DNS, and 
-0.4 to 0.35 for the TW, with a step of 0.15. The arrows indicate the cross stream velocity, 
scaled on magnitude (maximum 0.15(7). The middle row shows the streak structure over the 
wavelength of the TW, with contours of axial velocity at ±0.3 U (light/dark). The bottom row 
shows the axial vorticity, with contours at ±0.617/73 (light/dark). 
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Figure 9. hot and I uv (inset) as a function of time for TWs 2fe_0.625 (blue dashed line) and 
2a_1.25 (red solid line) starting at TW 46_3.125. Dots label the times of likely closest visits to 
each TW (times '2' and '3' coincide with those in figure^. 



exhibited turbulent episodes. Ideally, these should be composed from one very long - e.g. 
0(10,000 D/U) - turbulent run but, as discussed already, none could be generated at 
this Re where the TWs are available. As a result, two strategies at quantifying the TW 
visit frequency were undertaken. The first involved piecing together all the turbulent 
episodes produced during the suite of 5D runs described above. The second was to gen- 
erate longer turbulent data sets in a double-length 10D pipe; for relatively short pipes, 
the length of pipe appears to have a significant effect of the time turbulence is sustained 
for, as can be seen from Table These latter runs were almost exclusively initiated by 
randomly-selected velocity fields taken from a long turbulent coarse-grid run: see Table 
0] The coarse grid is too underresolved to reproduce the TWs accurately and hence cal- 
culate correlation functions directly but was perfectly adequate to initiate turbulence in 
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Figure 10. The surplus kinetic energy per unit mass in units of U 2 versus wall shear stress r 
m units of 2pU 2 /Re for udjvs starting at the upper branch TW 3&_3.125(+) (marked as a □ 
with a + in it). The laminar state is represented by the point (—4,0). All the TWs present are 
also plotted: A for 2-fold TWs, + for 3-fold TWs and o for 4-fold TWs. The numbered dots 
indicate the temporal points of closest approach. In chronological order: 1 - 36_2.5 (rightmost 
red circled +); 2 - 3j_2.5 (leftmost red circled +); 3 - 3e_2.5 (middle red circled +); 4 and 6 - 
2a_1.25 (large blue triangle); 5 and 7 - 26_1.875 (blue solid triangle). 



finer grid runs. 

For the 5D runs started at a TW, it was necessary to decide when the flow had left 
the neighbourhood of the TW and had started to behave in a turbulent fashion. Two 
different criteria were adopted for this. In the first (data set A), the start of the turbulent 
phase was taken as the time at which the correlation with the starting TW had reached 
a minimum. This turned out to be more stringent a criterion (i.e. produced a later time) 
than requiring that the flow trajectory merely enter the part of the KE — t plane occupied 
by the turbulent state. In the second (data set B), the start of the turbulent phase was 
taken as when the flow broke out of the initial rotational symmetry class of the starting 
TW (e.g. for run 4 this occurs at about 130 D/U: see FigureEJ). In both cases the end of 
the turbulent phase was well estimated by when the flow trajectory left this 'turbulent' 
K E — t area. The subsequent relaminarisation phase was easily identified with both the 
perturbation kinetic energy and wall stress decaying monotonically down to their lami- 
nar values (e.g. see figures l4l and ITU)l . These strategies produced a correlation data sets 
lasting 4, 154 D/U (A) and 2,363 D/U (B). For each, a joint probability density function 
was then computed for I to t and I uv corresponding to the TW with largest I to t at a given 
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Figure 11. I to t and I uv (inset) as a function of time for TWs 3&J2.5 (red solid line), 3c_2.5 
(black dotted line) and 3j_2.5 (blue dashed line) starting at TW 36_3.125. The dots label the 
time of closest visits (numbers correspond to points in figure 1X1)1 , 

time: see figures El an( l El There is a clear positive correlation between I to t and I uv 
and significant evidence for recurrent TW visits. Comparing the two pdfs, there are a 
lot of close visits during the interval after the initial transient but before the rotational 
symmetry class of the flow changes. This is presumably the result of the flow percolating 
out of the specific TW region in which it is initially inserted. As a result, the pdf from 
data set A is likely to overestimate the visit frequency. 

Using the criterion that a visit occurs if I to t > A & I to t + Iuv > 2A, the percentage 
visit time is plotted against the quality of the visit A in Figure El To assess how close 
a visit needs to be - i.e. what A is required - to be able to predict the instantaneous 
wall shear stress, Figures El and [21 plot the instantaneous DNS wall stress (over the 
matching wavelength) against the wall shear stress associated with the visited TW for 
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Figure 12. Comparison plots of the starting TW 3&.3.125 (top), the DNS flow (middle) and 
the TW 36_2.5 (bottom) at point 1 in figure tTTI The left column shows the velocity field at a slice 
where the DNS flow most matches 3&J2.5. The shading represents the axial velocity perturbation 
from laminar flow, with contours from -0.55 (dark) to 0.65 (light) (top), -0.55 to 0.5 (middle), 
and -0.55 to 0.35 (bottom), with a step of 0.15. The arrows indicate the cross stream velocity, 
scaled on magnitude (maximum 0.135 U). The right column shows the axial structure over a 
wavelength. Two contours of the axial velocity are shown at ±0.3 U (light/dark). 
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three different values of A. The choice of seeking maximal I to t to identify the closest 
TW was taken to facilitate the comparison as this should provide the best match be- 
tween the DNS streak structure at the pipe wall with that of the TW. Even with this, 
the correlation results for set A indicate that only A > 0.7 level visits are really good 
enough to start making predictions of the wall shear stress from the visited TW. Since 
there are no A = 0.7 visits in set B, A = 0.6 is used to indicate closer visits and better 
stress correlation is evident at this level compared to the A = 0.5 results. Since the ab- 
scissa is discretized over the 37 ttw values, each vertical strip of data in these figures 
indicates that that TW has been visited. Interestingly, Figure I2TJI shows a visit bias to 
TWs with larger wall stress (< —6 * 2pU 2 / Re) whereas Figure |2"T1 is oppositely skewed 
(> —6 * 2pU 2 1 Re) to lower wall stress TWs. A possible explanation for this is set A is 
dominated in a 2-to-l ratio with runs started by upper branch TWs compared to lower 
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Figure 14. The left plot shows the correlations Itot (upper line) and I uv (lower line) over one 
wavelength of TW 2a_1.25 at point 4 in figure lT^l The right plot shows the azimuthal distribution 
of the wall shear stress in units of 2pU 2 /Re at the axial position of maximum Itot and I uv near 
z = 1.5. The upper lines correspond to the azimuthal stress and the lower lines to axial stress 
(minus the laminar value of —4). The regular lines with 2-fold symmetry are for the TW and 
the more irregular ones from the DNS solution. 

branch TWs. In the initial adjustment phase where the flow trajectory gradually leaves 
the vicinity of the initial TW, the flow visits other TWs and these are more likely to 
be upper rather than lower branch solutions in the neighbourhood of an initial upper 
branch solution. The reason lower branch TWs feature more in set B may be because the 
turbulent episodes considered were not long enough to desensitise the visit statistics from 
the final relaminarisation process where the flow preferentially passes by lower branch 
TWs. This could indicate that efforts to remove this phase from the data but may not 
have been wholly successful. 

Further runs were carried out in a 27r/0.625 D (k, 10 D) pipe in the search for more 
sustained turbulent data. Extending the spatial domain by a factor of 2 allows a whole 
new set of TW wavelengths to fit into the pipe, roughly doubling the number of allowed 
TWs. The correlation calculations were not extended to encompass this enlarged set 
for two reasons. Firstly the new TWs are interspersed within the 5D set which already 
sample the available solutions well. Therefore the new TWs should have structures very 
similar to the existing TW set. Secondly, the increase in the numerical overhead: cur- 
rently calculating the correlation functions takes around 50% of the cpu time. However 
to compensate for this omission, the visit frequencies based on half the admissible TWs 
could justifiably be doubled. Initial conditions were randomly selected from an apparently 
sustained coarse turbulent run and used in intermediate- and fine-grid runs: see Table^J 
The statistics across the accumulated 5, 865 D/U long intermediate grid turbulent data 
were similar to those across the accumulated 4, 780 D/U long fine grid turbulent data so 
the sets were merged to give a one large data set over 10, 000 D/U in duration. The joint 
pdf of (Itot, Iuv) shown in Figure l2"4l resembles that from 5D data set B as does the wall 
stress comparison during visits shown in Figure [231 Again, the most visited TWs look 
to be those with lower wall stresses and essentially the same subset of TWs are visited 
closely - 2a_1.25, 26.1.875 and 36_2.5 in 5D data compared to 2a_1.25, 26.1.875, 36_2.5 
and 3c_3.125 in the 10D data. Figures l2"2l and l2"3l show an example of a visit to 2a_1.25 
during a typical turbulent episode in these suite of runs. Significantly, no 4- fold symmetric 
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Figure 15. Comparison plots of the DNS flow (left column) and the TW 2a_1.25 (right col- 
umn) at point 4 in figure EU The top row shows the velocity fields at the streamwise position of 
maximum hot shown in figure ITH f z ~ 1.5). The shading represents the axial velocity perturba- 
tion from laminar flow with contours from -0.55 (dark) to 0.5 (light), with a step of 0.15. The 
arrows indicate the cross stream velocity, scaled on magnitude (maximum 0.098 U). The middle 
row shows the streak structure over the wavelength of the TW, with contours of axial velocity 
at ±0.3(7 (light/dark). The bottom row shows the axial vorticity, with contours at ±0.8 [//D 
diaht/dark). 
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Figure 16. Comparison of the flow at point 5 in Figure HOI with TW's 26J..875 (left) and 
2a_1.25 (right). The top row shows the velocity fields for the TW's at the streamwise position of 
maximum correlation, the middle row the DNS flow, and the bottom row the wall shear stress 
at the same position. The shading represents the axial velocity perturbation from laminar flow 
with contours from -0.55 (dark) to 0.5 (light) (right column), -0.4 to 0.2 (top left), and -0.7 
to 0.65 (middle left), with a step of 0.15. The arrows indicate the cross stream velocity, scaled 
on magnitude (maximum 0.162 U). For the wall shear stress, the top lines show the azimuthal 
stress (in units of 2pU 2 / Re) , and the bottom lines axial stress (minus the laminar value of —4). 
The regular lines with 2-fold symmetry are for the TW's and the more irregular lines are from 
the DNS solution. 
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Figure 17. Joint pdf of Itot and I uv calculated over the turbulent episodes of all 5D runs started 
at the 8 upper and lower branch TWs (the values of I to t and I uv at a given time are selected 
by finding the TW with largest value of Itot) Contours are drawn at 0.01,0.1,0.3,0.5,0.7 and 
0.9 of the maximum value. The dashed lines cordon off the visit region defined by Itot > 0.5 & 

Itot + Iuv > 1. 

TWs were visited in the 10D runs reinforcing the observation made in the 5D work that 
these TWs are not in the same part of phase space as that populated by a turbulent flow. 

Figure ^5] brings together the results of the frequency analysis by comparing the visit 
percentage as a function of visit quality (A) for the three data sets. Doubling the 101? 
frequency data (as only half the possible TWs are considered) brings the observed visit 
frequency more into line with the results from the 51? data set B at A = 0.5. Taken 
together, they suggest that TWs are visited for approximately 10% of the time in tur- 
bulent pipe flow. Using Itot > A as the criterion for a visit would increase the frequency 
of visits, by, in effect, ignoring the cross-stream component of the flow and considering 
the streamwise (streak) structure only. 

Figure [2H shows the perturbation kinetic energy and the mean wall shear stress for a 
lOi? pipe for one of the intermediate and fine grid runs started from coarse grid turbulent 
data. This shows that the turbulent flow in this pipe is concentrated in a small patch of 
the KE — t space, with both grids occupying the same region of the space. The other 
intermediate and fine grid turbulent runs arc also concentrated in this region. Coarse grid 
turbulent runs occupy a somewhat different region of the KE — r space with lower values 
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Figure 18. Joint pdf of Itot and I U v calculated over the turbulent episodes of all 57? runs 
started at the 8 upper and lower branch TWs using time episodes after the flow leaves the 
initial symmetry class (the values of Itot and I uv at a given time are selected by finding the TW 
with largest value of I to t) Contours are drawn at 0.01, 0.1, 0.3, 0.5, 0.7 and 0.9 of the maximum 
value. The dashed lines cordon off the visit region defined by Itot > 0.5 & Itot + Iuv > 1. 



of both variables. Both the intermediate and fine grid have the same radial resolution 
(TV = 50). A test run was performed with the same grid in 9 and z as the intermediate 
grid but an increased resolution in s (N = 75). Again, the flow occupied the same region 
in the KE — r space. 

Figure |23 shows power spectra of the surplus kinetic energy and the wall shear stress 
for turbulent flow in a 10D pipe. The spectra were generated by taking the data for the 
four longest fine grid runs that were started from coarse grid turbulent data, with a sam- 
ple time of 782.1 D/U. A Hanning window was applied, with the windowed data scaled 
so that the variance matches that of the original data, and the spectra were generated by 
averaging over the four samples. Interestingly, the spectra show a peak, corresponding 
to a period of approximately 98 D/U which suggests a periodic orbit embedded in the 
turbulent dynamics f. This is currently under investigation. 



f The intermediate grid data produced spectra consistent with those for the fine grid, while 
the coarse grid also produced a peak, but at a higher frequency. 
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Figure 19. The percentage of time a turbulent flow 'visits' a TW based on the criterion hot > A 
and Itot + Iuv > 2A as a function of A. Statistics gathered from the 5D runs started at a TW with 
a) only the initial transient subtracted (dashed red line) and b) considering only times where 
the flow has left the initial symmetry class of the TW (solid red line with dots on), together 
with turbulent data compiled from various runs at WD (blue solid line with crosses on). 

4. Discussion 

At this point it is worth cataloguing the achievements of this investigation. 

(a) Using the travelling waves already known as starting points, all travelling wave 
solutions (TWs) of varying wavelengths and azimuthal symmetries which exist at Re = 
2400 have been traced out (Figures [T] and They naturally partition into 3 distinct 
classes of particular rotational symmetry about the axis - 2-fold, 3-fold and 4-fold sym- 
metric TWs - and 37 in total fit into a 7r/0.625 ( « 5) D periodic pipe. 

(&) The linear stability of four 'lower branch' TWs (TWs with some of the lowest wall 
stresses of the 37) and four 'upper branch' solutions (TWs with the some of the highest 
wall stresses of the 37) has been carried out. All are inertially unstable and therefore 
saddles in phase space with growth rates typically of the size 0(0.1 U/D) and all have 
very low dimensional unstable manifolds. 

(c) All the lower branch solutions considered appear to sit on a surface which separates 
initial conditions which uneventfully rclaminarise and those which lead to a turbulent- 
looking evolution. In contrast, initial conditions near all the upper branch solutions tested 



Recurrence of Travelling Waves 



27 



n r 



-10 



-11 - 



-12 



-®- 



i ! ! 



-9.5 -9 -8.5 -8 -7.5 -7 -6.5 -6 -5.5 -5 -4.5 -4 

x 

TW 

Figure 20. A plot of the wall shear stress tdns versus the wall shear stress ttw of the TW 
being visited (tdns is measured across the matching wavelength of the TW) for the 5D data 
set A runs. The visit criterion is I to t > A & hot + Iuv > 2A, where A = 0.5 (red dots) indicate 
approximate visits, A = 0.7 (blue circles) close visits and A = 0.9 (black squares) very close visits. 
The thick (black) diagonal dashed line indicates a perfect match tdns — ttw- The horizontal 
(blue) dashed line is the mean wall shear stress across the whole pipe, the dotted lines indicate 
one standard deviation either side of this mean and the limits of the vertical axis have been 
set to the maximum (-12.8) and minimum (-5.2) wall shear stress values (in units of 2pU 2 /Re). 
Each vertical strip of points represents one TW - the close visits only occur for upper branch 
solutions where ttw < —5.2. 



invariably become turbulent. 

(d) Turbulence in a 5D periodic pipe at Re = 2400 may be long-lived but ultimately 
appears only transient. Turbulence also seems transient in a 10D pipe at Re = 2400 but 
the transients have, on average, a longer life. The mean lifetime of a turbulent episode 
is sensitive to the numerical resolution used, with turbulence appearing to last longer, if 
not sustained, in undcrresolved calculations. 

(e) A number of different correlation functions were experimented with to measure 
how 'close' a given velocity field is to a TW. Two, I to t and were chosen and a visit 
criterion - I tot > A and Itot + Iuv > 2A - developed based on a 'quality of visit' param- 
eter A. After examining the velocity matches at various different levels of A, a value of 
0.5 was taken as indicating a 'visit' with values of 0.6 and higher indicating a 'close' visit. 
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Figure 21. A plot of the wall shear stress tons versus the wall shear stress ttw of the TW 
being visited {tons is measured across the matching wavelength of the TW) for the 5D data 
set B runs. The visit criterion is Itot > A & I to t + luv > 2A, where A = 0.5 (red dots) indicate 
approximate visits, A = 0.6 (blue circles) close visits. There are no very close (A = 0.9) visits. 
The thick (black) diagonal dashed line indicates a perfect match tons — ttw- The horizontal 
(blue) dashed line is the mean wall shear stress across the whole pipe, the dotted lines indicate 
one standard deviation either side of this mean and the limits of the vertical axis have been 
set to the maximum (-9.6) and minimum (-5.3) wall shear stress values (in units of 2pU 2 / Re). 
Each vertical strip of points represents one TW. For A = 0.5, the TWs visited are 2a_1.25, 
2&.1.875, 36_2.5, 3o_3.125 and 3c_3.125. The TWs closely visited (A = 0.6) are 2a_1.25, 2&.1.875 
and 36_2.5. 



(/) Turbulent data in both 5D and 10D pipes indicate that some 2-fold and 3-fold 
symmetric TWs are recurrently visited while others are not. In particular, no evidence 
was found for visits to 4-fold symmetric TWs. The visited TWs correspond to low-to- 
intermediate wall stress solutions which are embedded in the same part of phase space 
as the turbulent dynamics. Other TWs are clearly in very different phase space locations 
and hence are never visited except if the flow is specifically inserted there initially (e.g. 
all the 4-fold TWs). 

(g) Based on the correlation functions, Itot and I uv , used and the visiting criterion 
adopted (A = 0.5), numerical evidence suggests that travelling waves are only visited for 
about 10% of the time in turbulent pipe flow. 



The last finding answers the original motivating question for this study. The fact that 
a turbulent flow is apparently visiting a TW only 10% of the time implies that it is 
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Figure 22. The left plot shows the correlations Itot (solid line) and I uv (dashed line) over one 
wavelength of TW 2a_1.25 for a fine grid 10D pipe staring from a coarse grid turbulent flow. 
The right plot shows the azimuthal distribution of the wall shear stress in units of 2pU 2 /Re 
at the axial position of maximum Itot + Iuv near z — 1.84. The upper lines correspond to the 
azimuthal stress and the lower lines to axial stress (minus the laminar value of —4). The regular 
lines with 2-fold symmetry are for the TW and the more irregular ones from the DNS values. 

of limited use to view turbulence as the random switching between the neighbourhoods 
of TWs. From the perspective of, say, predicting the average turbulent wall stress, an 
appropriately weighted sum of all the relevant TW wall stresses seems unlikely to work 
given that so much time is spent away from the TWs in phase space. The fact that even 
during a visit, the match between the flows wall stress and that of the TW can be poor 
is further cause for pessimism. The way out of this conclusion, of course, is that other 
phase space objects such as periodic orbits - perhaps glimpsed in figure l?71 - need to be 
included in any such expansion. The results presented here indicate just how important 
these missing objects are given the low visit times for the TWs and the next challenge 
to isolate them seems clear. 
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Figure 23. Comparison plots of the DNS flow (left column) and the TW 2a_1.25 (right column) 
captured during a turbulent fine grid run in a WD pipe. The top row shows the velocity fields 
at the streamwise position of maximum hot + Iuv shown in figure 12^1 (z ~ 1.84). The shading 
represents the axial velocity perturbation from laminar flow with contours from -0.7 (dark) to 
0.35 (light) for the DNS, and -0.55 to 0.5 for the TW, with a step of 0.15. The arrows indicate the 
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Figure 25. A plot of the wall shear stress tdns versus the wall shear stress ttw of the 
TW being visited (tdns is measured across the matching wavelength of the TW) for the 10D 
runs. The visit criterion is Itot > A & I to t + Iuv > 2A, where A = 0.5 (red dots) indicates 
approximate visits and A = 0.6 (blue circles) closer visits (there are no A = 0.7 visits). The 
thick (black) diagonal dashed line indicates a perfect match tdns = ttw- The horizontal (blue) 
dashed line is the mean wall shear stress across the whole pipe, the dotted lines indicate one 
standard deviation either side of this mean and the limits of the vertical axis have been set 
to the maximum (-10.48) and minimum (-4.96) wall shear stress values (in units of 2pU 2 / Re). 
Each vertical strip of points represents one TW. For A = 0.5, the TWs visited are 2a_1.25, 
26.1.25, 3c_1.25, 2o_1.875, 26.1.875, 3c_1.875, 36_2.5, 3c_2.5, 3d_2.5, 3a_3.125 and 3c_3.125. The 
TWs closely visited (A = 0.6) are 2a_1.25, 26.1.875, 36_2.5 and 3c.3.125. 
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Figure 27. The power spectra of surplus kinetic energy (multiplied by 5000) in units of U 3 D 
(dashed) and wall shear stress in units of AUD/Re 2 (solid) versus frequency. The maximum in 
the spectra corresponds to a period of approximately 98 in units of D/U. 
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Table 2. Summary of visits for the 12 transitional runs (4 starting at lower branch TWs + 8 
starting at upper branch TWs) which have turbulent episodes (data set A). TWs visited using 
the criterion I tot > 0.5 & I to t + 4. > 1 at least once in a given run have the symbol o entered 
under that column against them (* indicates the starting TW). A • indicates where the higher 
threshold of I tot > 0.6 & hot + Iuv > 1.2 is satisfied. The initial transient as the flow moves away 
from the starting TW - typically of O(30 D/U) duration - is not considered. The code for the 
run numbers is: 1 - 26_1.25(— ) (where (— ) indicates the sign of the eigenfunction perturbation), 
2 - 3a_2.5(+), 3 - 37i_2.5(+), 4 - 46_3.125(+), 5 & 6 - 2o_L.25(+/-), 7 & 8 - 36_3.125(+/-), 9 
& 10 - 3j_2.5(+/-), 11 & 12 - 4c_3.125(+/-). 
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Table 3. As for Table except now only times when the flow has left the symmetry class of 
the starting TW are considered (data set B). Also numerical values are recorded of Itot + Iuv in 
all near-miss instances where hot > 0.5 but a o (A = 0.5) or • (A = 0.6) is not warranted. 
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Pipe Length Run Turbulent Duration Initiation Resolution 
in D/U (N, M, K) 
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433 (673) 


2b_1.25(-) 


(50,60,60) 
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327 (554) 


3a_2.5(+) 


(50,60,60) 
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340 (517) 


3/i_2.5(+) 


(50,60,60) 
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147 (251) 


4&_3.125(+) 


(50,60,60) 
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195 (385) 


2a_1.25(+) 


(50,60,60) 


6 


397 (609) 


2a_1.25(-) 


(50,60,60) 
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153 (265) 


36_3.125(+) 


(50,60,60) 
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263 (400) 


3fe_3.125(-) 
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(50,60,120) 
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(50,60,120) 


13 


558 
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(50,60,120) 


14 


265 
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(50,60,120) 


15 


640 
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(50,60,120) 


10 


1135 
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(50,60,120) 



Table 4. This Table lists the various runs used to produce the 'visit' statistics. The 5D data 
comes from initiating the code with a disturbed TW and two 'turbulent' duration times are 
listed. The first is measured from when the flow leaves the initial symmetry class of the TW 
(used for data set B) whereas the second, larger figure (in parentheses) is based on when the 
correlation function for the starting TW stops decreasing which is a much weaker condition 
that occurs earlier (used for data set A). The 10D runs are almost exclusively started using 
randomly selected velocity fields from a long turbulent run generated using a coarse (25,32,60) 
grid (a strategy labelled as 'T' in the Table). Note that some of the 10D runs were still turbulent 
when this data was gathered, whereas all the 5D runs had relaminarised. 



